Free planar isotropic-nematic interfaces in binary hard-rod fluids 
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Within the Onsager theory we study free planar isotropic-nematic interfaces in binary mixtures 
of hard rods. For sufBciently different particle shapes the bulk phase diagrams of these mixtures 
exhibit a triple point, where an isotropic (I) phase coexists with two nematic phases (A''i and of 
different composition. For all explored mixtures we find that upon approach of the triple point the 
IN2 interface shows complete wetting by an intervening A'^i film. We compute the surface tensions of 
isotropic-nematic interfaces, and find a remarkable increase with fractionation, similar to the effect 
in polydisperse hard-rod fluids. 
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I. INTRODUCTION 

Colloidal suspensions of rodlike particles are well- 
known to exhibit a wealthy phase behavior as a function 
of density 0, 0, Early experiments on vanadiumpen- 
toxide solutions jj| and suspensions of Tobacco Mosaic 
Virus (TMV) particles showed a first-order transition 
from a disordered (isotropic) fluid phase to an ordered 
(liquid-crystalline nematic) fluid phase upon increasing 
the concentration of rods sufficiently. The nematic is ho- 
mogeneous (it is a fluid), but the rods are on average 
oriented in a specific direction h, the so-called nematic 
director. The bulk nematic ordering is of uniaxial sym- 
metry, i.e., there is azimuthal symmetry about n. The 
isotropic-nematic (IN) transition was first explained by 
Onsager, who modelled the colloidal rods as hard nee- 
dles. He first showed that the average pairwise excluded 
volume is reduced in the nematic phase compared to the 
isotropic phase, and then argued that the resulting gain 
of free volume (and hence translational entropy) compen- 
sates the loss of orientation entropy (due to the nematic 
ordering) at sufficiently high concentrations of rods 0. 
In other words, the ordering follows as a consequence of 
maximizing the total entropy. More recently it was shown 
that a further increase of the concentration of TMV can 
also give rise to smectic ordering 0,13, where the transla- 
tional invariance is broken and the system forms a layered 
structure. In the 1980's and 1990's computer simulations 
0, , and later density functional theories 0, 0, ^3 ' 
have shown that also the smectic phase can be explained 
by the hard-rod model for colloidal rods. One concludes, 
therefore, that the bulk phase behavior of these systems 
is well understood by now. 

Onsager's theory for hard rods has been extended to 
describe bulk mixtures of colloidal rods. For the case of 
binary mixtures of longer and shorter rods, it was found 
that the IN transition is accompanied by strong fraction- 
ation, such that the coexisting nematic phase contains a 
relatively large fraction of the longer rods 0, 0| . Later 
theoretical work on long-short mixtures also showed the 
possibility of nematic- nematic {N1N2) demixing (driven 
by a peculiar competition between orientation entropy 
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and ideal mixing entropy), and 
nematic (IN-\N-2) triple point i 

[Til HtL ITa IT^ of mixtures with a length ratio more ex- 
treme than about 1:3. Later also binary mixtures of thin 
and thick hard rods were considered. They were shown 
to have phase diagrams similar to those of long-short 
mixtures, i.e., with strong fractionation at IN coexis- 
tence and with N1N2 and IN1N2 coexistence for diam- 
eter ratios exceeding 1: 3.8 [11 [19, 20, 21], but with an 
additional possibility for isotropic- isotropic ihli) phase 
coexistence due to the depletion effect [2^, [U [22| if the 
diameter ratio is more extreme than about 1:8. Interest- 
ingly, thin-thick mixtures have recently been realized ex- 
perimentally by mixing "bare" fd virus particles (length 
1/zm) with ones that are "coated" with polyethylenegly- 
col (PEG) [2^. The diameter ratio of these systems can 
be tuned by varying the ionic strength of the solvent: due 
to an increasing salt concentration the effective diameter 
of the (charged) bare rods shrinks because of enhanced 
screening, whereas that of the PEG-coated rods is not (or 
hardly) affected because of the steric nature of PEG. Ex- 
ploiting this effect allowed for the experimental study of 
diameter ratios up to about 1:4.5, and INi, IN2, N1N2 as 
well as IN1N2 triple coexistence were actually observed 

M 

The present study is devoted to the planar interfaces 
that exist between the coexisting bulk phases in binary 
mixtures of colloidal rods. Our focus is on the calcula- 
tion of both thermodynamic and structural properties of 
these interfaces. The main thermodynamic quantity of 
interest is the surface tension 7, and the structural prop- 
erties we will investigate are the profiles of the density 
and the order parameters. It is known from the study of 
the IN interface of pure (one-component) suspensions of 
rods that 7 depends on the angle between the interface 
normal z and the director n of the nematic phase asymp- 
totically far from the interface 24, 25, 26] . The study of 
Ref. [24L I27I ] showed that 7 is minimal when n ± z, and 
on this basis (and on the basis of some of our own test 
calculations) we assume this to be the case for mixtures 
as well. It is also established by now [^ [53, [13] that 
(i) the density profile and the nematic order parameter 



2 



profile of the IN interface of the pure hard needle fluid 
change monotonically from their values in the isotropic 
bulk phase to those in the nematic phase, (ii) the inter- 
face thickness is of the order of the length L of the rods, 
and (iii) the interfacial biaxiality is small and nonmono- 
tonic. In this paper we will show that the density profiles 
in mixtures of rods are not always monotonic, and that 
the interface thickness not always of order L due to the 
formation of macroscopically thick wetting films close to 
the bulk triple points. Some of these findings have been 
reported briefly elsewhere [s^l- Moreover, we will show 
that the tension in mixtures of rods tends to be sub- 
stantially higher than that of the pure systems of their 
components. 

This paper is organized as follows. In section ^ we 
introduce the Onsager functional and the basic Euler- 
Lagrange equation. In section UTTl we solve this equation 
for bulk geometries, and we present a few typical bulk 
phase diagrams. In section ITVI we present our method to 
solve the Euler-Lagrange equation for interface geome- 
tries, and study INi, N1N2 and IN2 interfaces, the lat- 
ter in particular in the vicinity of the bulk IN1N2 triple 
point. A summary and some discussion of our results will 
be presented in section Ivl 



II. DENSITY FUNCTIONAL 

We consider a fluid of hard cylinders of two different 
species ct = 1, 2 of length L^j and diameter D^j in a macro- 
scopic volume V at temperature T and chemical poten- 
tials /ia-. The thermodynamic properties and the struc- 
ture of the system can be determined from the grand po- 
tential functional r2[{pcr}] of the one-particle distribution 
functions /Ocr(r, w), where r denotes the center-of-mass co- 
ordinate of the rod of species a and Cj the orientation of 
the long axis. The functional is such that (i) it 

is minimized, for given ({/io-}, V, T), by the equilibrium 
one-particle distributions ^^(rjcZ'), and (ii) the minimal 
value of the functional is the equilibrium grand potential 

n [all. 

Within the second virial approximation and in the ab- 
sence of external potentials, the functional il[{pcr}] can 
be written 111 as 



mp^}] = [ dqpaiq)[HpAq)LlD,]-l-P^,. 

~\Y1 j d-qdq' f„„>{q]q')p„{q)p„,{q'), (1) 



with (3 = (kT)^^ the inverse temperature and faa'{q',q') 
the Mayer function of the acr'-pair of rods with coordi- 
nates q = {r, (I)} and q' — {r',(I)'}. For hard rods, the 
focus of our study, faa'{q',q') equals —1 if the rods over- 
lap and vanishes otherwise. Onsager argued that the sec- 
ond virial approximation is accurate for long rods, and 
becomes even exact for isotropic and nematic bulk flu- 
ids in the limit of vanishing diameter-to- length ratio 



We shall adopt this limit throughout this paper, i.e., we 
consider D„/L„i for any aa' pair. Therefore, the 
relative shape of the rods is only characterized by the 
ratios I — L2/L1 and d — D2/ Di of the lengths and the 
diameters, respectively. 

The minimum conditions Sfl[{p^}]/ d pa{q) = on the 
functional lead to the set of nonlinear integral equations 

Hp^iq)LlD,] - E/ dq'Ua'{q;q')Pa'(q') = (2) 



to be solved for the equilibrium distributions p^ (q) . Once 
determined, they can be inserted into the functional to 
obtain the equilibrium value of the grand potential 



(3) 



P^^\Y.j dqpAq){HpAq)LlD,]-2 



Note that = —pV for a bulk system in a volume V , with 
p — p{{pa}, T) the pressure. In the presence of a planar 
surface or interface of area A we have VI — —pV + 7A 
with 7 = 7({/Xct},T) the surface or interface tension. 

In general, fluctuations of the interface position (cap- 
illary waves) are important in the analysis of fluid-fluid 
interfaces. The amplitude of such fluctuations is con- 
trolled by the "wetting parameter" [33.13^ 



47r7^2 



(4) 



where ^ is the bulk correlation length. For rods of typi- 
cal length L and diameter D we shall see that f ~ L and 
7 ~ ksT/LD, and hence uj ~ D/L, i-e., uj vanishes in 
the Onsager limit. As a consequence the capillary- wave 
fluctuations are unimportant, i.e., the mean-field density 
functional is sufficient to describe interfacial phenom- 
ena in fluids of long hard rods. 



III. BULK PHASE DIAGRAMS 

The bulk thermodynamic properties of binary hard- 
rod fluids were studied extensively within Onsager the- 
ory. The minimization of the functional was either per- 
formed variationally [ll Elfiljilpi or through 
a fully numerical solution |14L I21L l28l | . We adopted the 
latter approach as it can be easily generalized for inho- 
mogeneous systems. For clarity we briefly repeat the es- 
sential points of the method and summarize the available 
results. 

The bulk distribution functions of the isotropic and ne- 
matic phase are translationally invariant, i.e., Pa{Y,ijj) = 
Pcr{(jJ), which allows to reduce Eq. ((2Jl to 



HpA^)LlDa 



E 



(5) 
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with E^^i the excluded volume of a pair of cylinders of 
species a and a' given by 

= L„L,,{D, + D,,)\smj\ (6) 

in terms of the angle 7 between (2; and a)', i.e., 7 = 
arccos(a) • a)'). Note that additional 0{LD^) terms are 
being ignored in Eq. JSJ, in line with the needle limit 
{Da /La' — > 0) of interest here. 

At sufficiently low {Pficr} the only stable solution of 
Eq. (0 is the isotropic distribution = Ua/i'in), 

with Ha- = J dwpa-{tij) the bulk number density of species 
a. As /Iter are high enough, one or, possibly, two sets 
of stable uniaxial solutions {(^) = Pa{0) exist, with 
9 = arccos(a) • n) the angle between uj and the nematic 
director -h. These distributions have "up-down" symme- 
try. Pa {6) — Pa{T^ — 0), hence Pcr(0) needs only to be 
determined for 9 £ [0,7r/2]. Using an equidistant 6'-grid 
of Ne = 30 points 9, e [0,7r/2], where 1 < i < iV^, we 
iteratively solve Eq. for the set of 2Ne equations in or- 
der to find Pa{9i) numerically. The integral in Eq. (|SJ) is 
calculated with the trapezoidal rule. Coexistence of dif- 
ferent phases {/, A'^i, N2} can be determined by imposing 
conditions of mechanical and chemical equilibrium. 

In order to gauge the accuracy of the chosen 6'-grid we 
calculate the resulting densities of the coexisting isotropic 
and nematic phase of the one-component system. We find 
n'L^D{-K/4) = 3.281 ± 0.001, n^L^Din/A) = 4.172 ± 
0.001; the nematic order parameters of the two coexisting 
phases are = 0.008 ± 0.001 and = 0.791 ± 0.001, 
and the pressure {■k/4)PpL'^D = 14.045 ± 0.001. These 
data, based on N0 = 30, differ by less then a percent 
from the most accurate results available in the literature 
0, which we can reproduce with Ng > 80. In order to 
have consistency between bulk and interfacial results we 
take Ng = 30 in most of our calculations. The exception 
is the case of long-short mixtures in the nematic phase, 
which requires Ng — 50 for acceptable accuracy. 

In Fig. ^we show both pressure-composition (a) and 
density-density (b) representations of bulk phase dia- 
grams of thin-thick binary mixtures {La = L, D2 > -Di) 
for several diameter ratios d = D2/D1. In Fig. ^a) the 
composition variable x = ^2/(711 + 712) denotes the mole 
fraction of thick rods and / = {p-pthick)/ {pthin -Pthick) 
is a dimcnsionless shifted pressure which takes the values 
f — 1, at isotropic-nematic coexistence of the pure-thin 
{x = 0) and pure-thick {x = 1) systems, respectively. 
Note that {7^/4)ppth^nL''Dl = {TT/4)PptMckL^D2 = 
14.045, i.e., pthick = Pthin/ d, and that the tie-lines con- 
necting coexisting phases are horizontal in the f — x rep- 
resentation of Fig. da). 

At low pressures (or low densities) the phase diagrams 
show an isotropic (/) phase and at higher pressures (or 
densities) one {d = 3.0, 3.5) or two {d = 4.0, 4.2) nematic 
phases (A^i and N2)- For diameter ratios d = 3.0, 3.5 the 
phase diagram is spindlelike, and the only feature is a 
strong fractionation at coexistence, such that the nematic 



phase is relatively rich in thick rods and the isotropic 
phase in thin ones. The reason behind this fractionation 
is the relatively large excluded volume in interactions 
of the thick rods, which makes them more susceptible 
to oricntational ordering 0, IT^ |21.] . The fractionation 
of isotropic-nematic coexistence becomes stronger for in- 
creasing d. For 3.8 < d < 4.29 the bulk phase diagram 
develops nematic-nematic {N1N2) coexistence in a pres- 
sure regime pt < p < Pc, with pt the triple-point pressure 
and Pc the critical (consolute) pressure. For d = 4.0 the 
consolute point is indicated by (*) in Fig. ^ The mecha- 
nism of this demixing transition was spelled out in detail 
m Refs.[ll|23, and involves a competition between ori- 
cntational entropy (favoring demixing) and entropy of 
mixing (favoring mixing). Interestingly, the width of 
fractionation gap Ax = xmi — xj^i^ for the triple point 
A^i and N2 phases scales linearly with the triple-point 
pressure pt- The critical pressure of the N1N2 transition 
diverges as d ^ 4.29 J.7, 2lj. For d > 3.8 the lower 
bound of the N1N2 coexistence is an IN1N2 triple point, 
indicated by triangles (A) in Fig. ^for d = 4.0. With 
increasing d the triple point / and iVi phases approach 
the pure-thin bulk coexistence (i.e., xi,Ni — > 0), whereas 
the composition of the triple point 7V2 phase shifts to a 
pure-thick phase (x^Va ^ 1)- 

The f — X representation is convenient for our analy- 
sis, whereas the densities (volume fractions) of thin and 
thick rods are experimental control parameters For 
this reason the same phase diagrams of thin-thick binary 
mixtures are shown in Fig. in density-density repre- 
sentation, with n* = niLD'l{TT / A) and = niL£)|(7r/4) 
being the dimcnsionless bulk number densities of thin 
and thick rods, respectively. In this representation the 
tie-lines are no longer horizontal. 

In Fig. 12 a set of bulk phase diagrams for long- 
short binary mixtures {Da = D,L2 > Li) for sev- 
eral length ratios I = L2/L1 is presented. Here x = 
n2/{ni + 7^2) denotes the fraction of long rods and / = 

(P - Plong)/{p short " Plong) with piong = Pshort/l"^- AH 

characteristic features of the phase diagrams are the same 
as in thin-thick mixtures. The fractionation of the coex- 
isting IN2 and N1N2 phases has a strong dependence 
on I, and limits the values accessible for calculations to 
I < 3.1 for the chosen grids and the required accuracy. 
The main reason is that the nematic ordering in the triple 
point A'^2 phase is very pronounced, requiring a fine grid 



IV. FREE INTERFACES 

We now turn to the thermodynamics and the struc- 
ture of the free interfaces between the coexisting phases. 
We assume that the interfaces are planar, with surface 
normal z. The nematic director n of the asymptotic ne- 
matic bulk phase(s) can, in general, have a nontrivial tilt 
angle 9t = arccos(n • z) with respect to the interface nor- 
mal. In the present calculations we restrict attention to 
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9t = 7r/2, i.e., n _L z. This geometry is known to be ther- 
modynamically favorable because of its minimal surface 
tension [2ll2a|. 

The equilibrium distribution functions Paiz, i^) depend 
on the spatial coordinate z = z ■ r, and angular co- 
ordinates uj — {0,ip) defined by cos 9 — h ■ lu and 
sin0sin(/3 = £•(£>. These functions are solutions of the 
Euler-Lagrange equations Eq. Q at the coexisting chem- 
ical potentials Pfi^ = f3f^%°'^^, with boundary conditions 
Pa-{z zLco^Lo) = p^\d) being the two coexisting bulk 
distributions (labelled by and (— ) here for brevity). 

The planar symmetry, i.e., the independence of 
of the in-plane coordinates x and y, allows for 
a reduction of the numerics, since the "excluded slab" 
ICaa'{z,u;,z',6j') = — J dx'dy' f„ni{Y,uj\Y' ,uj') can be cal- 
culated analytically |25l |28| . This reduces the Euler- 
Lagrange equations to 

(3pT'- = \n[p,{z,Q:)LlD,]+Y. j dz' J du' 

a' 

XJC„„, {\Z - z'l, Uj')pa> (z', Lj'), (7) 

where the expression for /Co-o-' (|z|, a), a)') is given in Ap- 
pendix. 

In principle one could now solve Eq. l(7|l on a 
(z, 6, (f) grid. However, the numerical efforts can be fur- 
ther reduced if one realizes that biaxiality, i.e., the (p- 
dependence, is weak j22i,,2^]. In that case the truncated 



expansion 

M 

Pa{z,e,ip) = ^ p„^rn{z,0)cOs{2m,(p) (8) 
m=0 

is expected to be accurate for small M, and hence only a 
few "coefficients" pcr,m{z,6) {m < M) need to be deter- 
mined on a {z^ 0)-grid. It is important to realize, however, 
that Eq. JHJ implicitly assumes that the nematic direc- 
tor n does not vary in space. The coefficients Pa;miz-,0) 
follow from an insertion of Eq. ((HJ into Eq. , multipli- 
cation by cos(2to(^), and integration over ip {0 < ip < 2ti). 
For M = this yields 

= \n[p,,Q{z,6)LlD,]+Yl j dz'dO' sinO' 

a' 

xlCZ{z-z\9,9')p,,,o{z',e'), 

(9) 

where /C™ , (z- z', 0, 0') = (2^)-^ J^"" J^"" dipdip'IC,^,{z- 
z',uj,uj') is the doubly azimuthally integrated excluded 
slab, which we determine numerically once on an appro- 
priate grid. 

The lowest-order correction that takes into account 
biaxiality results from M = 1, and yields pa.m{z,9), 
(to = 0, 1) from the coupled set of equations 



= lnM^,0)LlD,]+Io 



P<y,iiz, 



^Pa,o{z,0), 

j dz'de' sin 6' (k.Z' {z - z', 9, 9')p^,,o{z', 9') + K.^, (z - z', 9, 9')p„,,^{z\ 9')) , 



= /: 



P<y,l{z, 



^ j dz'd9' sin 9' (k.I°„, {z ~ z\ 9,9') Pa', 9') 



Pafi{z,9)_ 

+KZ,{z-z',9,9')p,,,,{z',9')), (10) 

with KJ^J^,{z — z\9,9') — {2n)^^ J^^ J^^ dipd(p' cos{2k(p) coa{2nnp')ICcrcr' {z — z',lu,lj'), k,m = {0,1}, again to be 
determined numerically only once, with 



1 Z"^'' 1 + Vl - 

Io{x) = 7j- dip\ii(l + xcos(2ip)) —In- 



27r 7o 2 



Ii{x) — — / (iiy9 cos(2(^) In (l + a; cos(2(^)) 
2n Jq 

for |x| < 1. 

Note that the boundary conditions imply that \pa,i{z,9)LaD^\ 1. In the remainder of this paper we 

Pa-,i{z,9) for |z| oo. In general, the solutions shall mainly concentrate on solutions of Eq. (Q, although 

P<7,o{z,9) for M = are not identical to pa-,o{z,9) for some results of Eqs. (|10ll will be discussed. 

M = 1, but the difference is small in most cases since gy iteration of Eq. © (or Eqs. CS) with the appro- 



5 



priate boundary conditions we calculated pafi{z,0) (and 
Pa,i{z, 9)) for a number of state points on the INi, 

IN2 and A''iA^2 binodals. We used an equidistant spa- 
tial grid of Nz — 200 points Zi € [— 5L, 5i], an equidis- 
tant angular grid of Ng = 30 points 8j E [0, 7r/2] for 
thin-thick mixtures or an angular grid of No = 50 points 
Oj G [0, 7r/2] for long-short mixtures. From the equilib- 
rium distributions po-.oC-z, &) we calculated the local den- 
sity and the nematic order parameter profiles 

na-{z) = Att d9smdpcrfl{z,9) 



Sa{z) = 47r / d9sv[\9P2{cos9)pafi{z,9)/na{z), 
Jo 

with ^2(2^) = (3a;^ — l)/2 the second Legendre polyno- 
mial. In the case of iterating Eqs. Hl()|l the biaxiality is 
defined as El El 



A,(z) 



O 9 

- sin 9 cos 2o9 
2 ^ 



3 r^^ 

— / d9 sinH p^^iiz, 



The interface thickness t is defined as t = |z+ — z_| 
where z± are solutions of ni"'{z) = 0, where a prime 
denotes differentiation with respect to z. As this equa- 
tion has a set of solutions in every interfacial region, we 
choose {z±) be the outermost ones, i.e., the nearest to 
the bulk phases. This criterion provides a single measure 
for the thickness of both monotonic and non-monotonic 
profiles, with and without a thick film in between the 
asymptotic bulk phases at z — s- ±00. Also thin (or short) 
rods have a smaller excluded volume and a non- vanishing 
concentration in both coexisting phases, so their den- 
sity is a convenient representation of structural changes 
within the interface. The interfacial width for the one- 
component IN interface is, with the present definition, 
given by t/L = 0.697. We have checked that other defi- 
nitions of the thickness lead to similar results. 



A. INi and N1N2 interfaces. 

The /A''i interfaces exist only in a small pressure regime 
Pthin ^ P > Pt- They closely resemble the IN interface 
of the pure hard-rod fluid, i.e., the profiles of the order 
parameters Sa-{z) and the densities na-{z) change mono- 
tonically from the bulk values in the / phase to those in 
the A^i phase. 

The thickness of the INi interface is of order L which is 
similar to that of the pure system. With increasing d the 
INi surface tension at triple phase coexistence {p — pt) 
decreases monotonically to the IN surface tension of the 
pure system, as shown in Fig. |21 where the dimcnsionless 
surface tension 7* = is plotted. This is expected 

from an inspection of the phase diagrams as x/.atj — > 



for increasing d. For the diameter ratio d = 4.0 the 
surface tension at p = p( is given by 7^^^ = 0.209 ± 0.001 
with the same level of accuracy for all other calculations 
of the surface tensions. For the long-short mixtures the 
behavior of the surface tension at the INi interface as a 
fimction of length ratio I is very similar. For I — 3.0 the 
surface tension at p = pt is given by 7|jv^ = 0.212±0.001. 

In general, the order parameter and density profiles 
are shifted with respect to each other. Such a shift can 
be characterized by the distance 5 = \zn — zs\ between 
the centers z„ and of the density and order parameter 
profiles, defined by |2J| 



n{zn) 
S{zs 



1 , 



where +/— indicate asymptotic bulk values. A non-zero 
shift 6 reflects the fact that the thickness of the inter- 
face is different for rods of different orientations. For 
monodisperse rods it was found that S — 0.50L |2j. For 
binary mixtures S can be determined for each component 
separately. For thin-thick mixtures with d = 4.0 the INi 
interface at triple phase coexistence shows Sthin — 0.35i 
and Sthick = 0.55L. For long-short mixtures with I — 3.0 
the effect is similar for the short rods {Sshort — 0.37Li) 
and much more pronounced for the long rods {6iong — 
1.54Li, i.e., (5;o„g ^ 0.51L2)- 

The profiles of Saiz) and na-{z) at N1N2 interfaces 
are also monotonic. For the diameter ratio d — 4.0 at 
p — Pt the interface thickness is given by t/L — 0.592 
and the surface tension 7^^^.^ = 0.019 ± 0.001, which 
is an order of magnitude smaller than 7jjyj^ . Upon the 
approach of the critical point t/L — > 00 and the surface 
tension vanishes. For d > 4.0 surface tension 'j'niN2 (^^ 
p = Pt) increases approximately linearly with d as shown 
in Fig. El 

The biaxiality is found to be small in both the INi 
and the N1N2 interfaces. In Fig. 21 we present the profiles 
Aa-{z) of the triple point INi and N1N2 (inset) interfaces 
of the thin-thick mixture with d = 4.0, as well as that 
of the IN2 interface to be discussed later. The marked 
curves represent A2(z) (thick rods), the unmarked ones 
Ai(z) (thin rods). Figure H reveals that \Aa{z)\ < 0.017 
in the INi interface, and |Acr(z)| < 4.0 • 10""' in the 
N1N2 interface. Such small biaxialities indicate that the 
expansion of Eq. jSl truncated at Af = 1, is accurate for 
calculating Ag- (z) , while a truncation at M = yields ac- 
curate tensions and density profiles. In fact, we checked 
that the difference between the tensions based on uni- 
axial (A/ — 0) and biaxial {M — 1) profiles falls within 
the numerical accuracy, i.e., less than 1%. Our numerical 
data for the INi interface is consistent with that of Rcf. 
[23,^^. Even though the magnitude of Acr{z) is small, 
it is interesting to consider the structure of the profiles in 
some more detail. The first observation we make is that 
Act(z) > (< 0) at the isotropic (nematic) side of the 
INi interface for both species a = 1,2. This indicates 
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that rods at the /-side of the interface have a (small) 
preference for splay in the XY-plane, whereas those at 
the Ni side tend to "stick" through the interface (into 
the I side). A similar effect exists in the N1N2 interface, 
but now both species have an opposite tendency (see in- 
set): the thin rods splay in the XY plane at the Ni side, 
whereas the thick ones "stick" through, and vice versa at 
the N2 side. Recall, however, that these effects are small. 



B. The IN2 interfaces. 

The IN2 interfaces exist, for d > 3.8, in a pressure 
regime Ptiuck < P < Pt- The properties of the IN2 in- 
terfaces depend strongly on the pressure difference with 
the triple-point {IN1N2 phase coexistence). The surface 
tension of the IN2 interface shows a non-monotonic de- 
pendence on the bulk pressure p. It develops a maximum, 
which is several times larger than a linear interpolation 
between the tension of the two pure systems, as shown 
in Fig. 13 It turns out that the non-monotonic character 
of "iiN-iiv) is related to the fractionation at the IN2 co- 
existence, i.e., a larger composition change through the 
interface leads to a larger interfacial tension. However, 
the surface tension which corresponds to the pressure of 
maximal bulk fractionation (indicated by the dashed line 
in Fig. |SJ) is lower than the maximum of Jin2{p)- The 
maximal interface stiffness grows with species diameter 
ratio d as the composition difference between / and N2 
phases increases (see Fig. ^ . We have also compared the 
maximal surface tensions for different orientations of the 
director (n _L z and n || S) in several thin-thick mixtures 
and found the geometry n _L i to be thermodynamically 
stable, i.e., jn±z < lfi\\z by at least a factor of two. 

The relatively large surface tension of a mixture of rods 
compared to that of the pure systems of its components 
may well be an explanation for the relatively large ten- 
sions that were measured in suspensions of cellulose [s^ , 
which are known to be very polydisperse. This remains 
to be investigated in detail, however. 

The dimensionless undersaturation e = 1 — p/pt is a 
convenient measure of the pressure difference with the 
triple point. For 0.01 < e < 0.5 the profiles of the order 
parameters S^{z) and the density of the thick component 
n2{z) are smooth and monotonic, whereas ni{z) shows an 
accumulation of thin rods at the isotropic side of the in- 
terface. This effect becomes more pronounced for small 
undersaturation, i.e., e — s- 0, when a film of the nematic 
phase A^i appears in the IN2 interface. Note that the A^i 
phase is a metastable bulk phase for any e > 0, so the 
film thickness is finite. For d = 4.0 several profiles ni{z) 
and n2{z) for different values of e are presented in Fig. 
|S1 which clearly shows the film formation when e ^ 0. 
The asymptotic densities at z — > ±00 are those of the 
coexisting I and N2 bulk phases (at the corresponding 
e). Using translational invariance of the interface be- 
tween the bulk phases, we have shifted the profiles with 
respect to each other such that their INi interfaces co- 



incide. This shows that the local density of thin (thick) 
rods in the growing film remains constant, and exactly 
corresponds to the thin (thick) -rod density of the bulk 
triple point phase A'^i (indicated by the dashed lines in 
Fig. (HJ. The same identification can be made for all 
d (or / for long-short mixtures) as well as for the order 
parameter profiles Sa-{z). 

The biaxiality of the IN2 interface was found to be 
small. A typical profile for thin-thick binary mixture 
with d = 4.0 at e = 5 • 10"'' is presented in Fig. 2| The 
IN2 biaxiality profile can be considered as a composition 
of the (earlier discussed) INi and N1N2 profiles which is 
expected as the thickness of the wetting A^i film is larger 
than L. 

For all explored mixtures the thickness of the interface 
t/L (or the adsorption) was found to diverge logarithmi- 
cally with e ^ as shown in Fig. [7| For short-ranged 
interactions one expects, on the basis of mean-field the- 
ory [s^l, that 



(11) 



where C is an irrelevant constant offset, and ^ is the 
correlation length of the wetting film. This implies that 
the bulk correlation length ^at^ of the wetting Ni phase 
should follow from the slope of the logarithmic growth 
of t/L in Fig. □ The asterisks (*) in Fig. IHlshow the 
resulting values of S^Ni as a function of d. These can be 
compared to the correlation length that one can extract 
from the asymptotic decay of the one-particle distribu- 
tions Pa{z, uj) of the INi interface into the bulk A^i phase. 
This decay can best be analyzed in terms of the deviation 
from the Ni bulk density Spa{z,uj) — pa{z,il!) — p^^{bj), 
which we find to decay as 



5pc,{z,Lb) = ylo.(a')exp(-z/^ArJ, 



00, 



(12) 



where the only species and orientation dependence is in 
the decay amplitude Ac^{u)), i.e., the decay (correlation) 
length <^ is one and the same for all species and orien- 
tations. Such a "decay-law", with a single correlation 
length, is well-known in mixtures of simple liquids 3(|. 
The form H12|l is illustrated for d = 4.0 at p ~ pt in the 
inset (a) of Fig. |S1 where all curves (representing different 
6''s) are parallel on a logarithmic scale. The correlation 
length follows from the (common) slope of these curves, 
and is marked by (o) in the main figure. The agreement 
with the values of ^/jvi obtained from the logarithmic 
growth of the interface thickness is clearly good. The 
inset (b) of Fig. |S| shows the similar dependence of the 
bulk correlation length of triple-point A^i phase in the 
case of long-short mixtures, i.e., as a function of rods 
length ratio I, using the same symbols. 

In order to verify the thermodynamic condition of com- 
plete wetting, ^iN2 — liNi + at the triple-point 
pressure p = pt, we determine the ratio of surface ten- 
sions 



7/JV2 (e) 



limpipt (7/Ari +IN1N2)'' 



(13) 
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shown in Fig. |51 For all diameter ratios d considered here 
lim£_,o -R(e) = 1, which implies a vanishing contact angle. 
This constitutes the thermodynamic proof of complete 
triple point wetting in all thin-thick hard-rod mixtures 
with d > 4. For mixtures of long-short rods the behavior 
of i?(e) is the same, as shown in the inset of Fig. 
A mean-field analysis of the asymptotic behavior of the 
surface tension (e — > 0) in the case of complete wetting 
shows that [s^ 

i?(e) - 1 ~ e"^'" (14) 

with the critical exponent a — 1. Analysis of our results 
in Fig. I^lgives a = 1.00±0.05 which can be considered as 
a rough consistency test of our mean-field calculations. 

V. SUMMARY AND DISCUSSION 

We have studied free interfaces of binary mixtures of 
hard rods of either different diameters or different lengths 
within Onsager's second virial functional. On the basis 
of a vanishingly small wetting parameter, which implies 
that capillary fluctuations are not important, we argued 
that this mean-field functional provides a realistic de- 
scription of isotropic-nematic interfaces of long hard rods. 
We focussed on diameter ratios d > 3.8 (and length ra- 
tios I > 3.0), for which the bulk phase diagram exhibits 
an IN1N2 triple point, and restricted attention to the 
case where n_Lz, with n the bulk nematic director and 
z the interface normal. This is the thermodynamically 
most favorable geometry. 

We have determined the behavior of the surface ten- 
sions of /-/Vi, N1N2 and IN2 interfaces between coex- 
isting isotropic and nematic phases as a function of the 
bulk pressure and/or the diameter ratio d (length ratio 
I, respectively). The tension jini is always very close to 
the tension of the pure fluid of thin (or short) rods, and 
"fNiN2 varies from zero at the consolute N1N2 point to 
values as large as Jini at the triple point for d ~ 4 — 5. 
The thickness of the INi and N1N2 interfaces are always 
of the order of the lengths of the rods, except close to 
the N1N2 consolute point, of course. The surface tension 
7/JV2 is found to change non-monotonically with pres- 
sure, exhibiting a maximum close to (but not at) that 
pressure where the bulk fractionation is maximal. This 
maximum surface tension is considerably larger than that 
of the pure systems of the components, typically by a 
factor of order 3-5, not unlike the findings of Ref. |28l| . 
where the case d — S (without any triple point) was 
studied. The biaxiality was found to be very small in 
all cases, similar to the findings in Refs. l28| for the 
one-component case. Perhaps our most interesting find- 



ing is the phenomenon of complete triple-point wetting 
of the IN2 interface by an A^i film. The thickness of 
this film is found to diverge as — ^ln(l — p/pt) when 
p ^ Pt, with ^ the correlation length of the bulk A^i 
phase, pt the triple-point pressure, and p < pt the pres- 
sure. The triple-point wetting phenomenon is confirmed 
by the numerical value of the surface tensions, which sat- 
isfy limp-fpj 7/JV2 {p) = jiNi (pt) + IN1N2 (Pt)- Such a com- 
plete wetting scenario was found for all diameter ratios 
3.9 < d < 5.2 and length ratios 2.9 < ^ < 3.1 studied 
here. We expect that this finding will also hold for more 
extreme ratios c? > 5.2 and I > 3.1, which are more dif- 
ficult to analyze numerically because of the pronounced 
nematic ordering of the triple-point N2 phase. 

The predicted phenomenon of triple-point wetting may 
well be observable in the experimental system of bare 
and PEG-coated fd virus particles |2^ mentioned in the 
Introduction. We hope that this work stimulates further 
experimental activities in this direction. 

Another interesting direction for future theoretical 
work would be to consider isotropic-nematic interfaces 
of polydisperse mixtures, e.g. extending the theory for 
bulk systems developed in Ref. (s^l- For suspensions of 
length-polydisperse cellulose experimental measurements 
of the surface tension have been performed 35], and show 
that the surface tension is much larger than that of a pure 
system of rods. It is tentative to speculate that the frac- 
tionation effect that is also present in binary systems may 
explain this increase of the tension. We hope to address 
this question in future work. 
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APPENDIX 

In terms oi A = ^max{LrjUJ ■ z,L^'6j' ■ z), B = 
i min(Lo-'^ • z, L„ilj' ■ z) and the excluded volume Ea^i = 
L„L„i {D„ + -Do") I sin(arccos(w ■ ^JJ'))\^ the results of Ref. 
[23 reduce for Da/La — > to the following expression 
for the "excluded slab" used in Eq. 
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r 0, 



A^aa'(k|,W,w') = < 



4\AB\ 

E{UJ,LJ') 

2\A\ ' 



\z\ > \A\ + \B\, 
{\A\ + \B\ - \z\), \A\ - \B\ < \z\ < \A\ + \B\ 

\z\<\{\A\-\B\)\. 
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FIGURE CAPTIONS 

1. (a) Bulk phase diagrams of binary thin-thick mix- 
tures for different diameter ratios d in the f — x rep- 
resentation, with / = (p-pthick) / {pthin-Pthick) the 
dimensionless shifted pressure, and x the mole frac- 
tion of the thicker rods. We distinguish the fully 
symmetric isotropic phase (/) and orientationally 
ordered nematic phases (iVi and A^2)- For the di- 
ameter ratio d = 4.0 the IN1N2 triple phase coexis- 
tence is marked by (A), and the N1N2 critical point 
by (*). (b) The same phase diagrams in density- 
density representation, where = niLD\{-K/A) 
and = n2LD2{'K / ^) are the dimensionless bulk 
number densities of thin and thick rods, respec- 
tively. The tie-lines connect coexisting state points. 

2. Bulk phase diagrams of binary long-short mixtures 
for different length ratio / in the f — x represen- 
tation (see caption to Fig. P). (A) mark IN1N2 
triple phase coexistence and (*) marks N1N2 criti- 
cal point for / = 3.0. 

3. Dimensionless surface tension 7* = (3j/LDi of INi 
(o) and N1N2 (o) interfaces at triple phase coexis- 
tence {p = Pt) for different diameter ratio d of thin- 
thick mixtures. The dashed line corresponds to the 
surface tension of the one-component IN interface 
for which 7|jv = 0.156 ± 0.001 

4. Biaxiality profiles ^a{z) of the thin (without sym- 
bols) and thick (marked by x ) rods in the INi and 
IN2 (at undersaturation e = 5 • lO^**) interfaces for 
diameter ratio d — 4.0. The inset shows the same 
quantity for the N1N2 interface. 

5. Dimensionless surface tension "f}j^^ = P"/in2/ LDi 
at IN2 interfaces as a function of dimensionless 
pressure p* = PpL'^Di{tt/A) for different diameter 
ratio d= 4.0(0), 4.2(o), 4.5(<i), 5.0([>) of thin-thick 
mixtures. The dashed line indicates the pressure of 
maximum fractionation. The data for d = 3.0(0) 
are included for comparison (from Ref. 28]). 

6. Density profiles of the thin rods nl{z) (a) and the 
thick rods 712(2) (b) in the IN2 interface for di- 
ameter ratio d = 4.0 at triple point undersatura- 
tions e = 1 - p/pt = 0.29, 0.1, 0.01, 5 x lO"'^, 1.3 x 
10"'',2.5 X 10-5. The bulk I/N2 phase is at 
z — )■ — cx)/cx3. The dashed lines nl — 3.977 and 
11,2 — 0.312 represent the bulk density of thin 
(thick) rods in the triple point A^i phase. These 
profiles indicate the formation of a wetting A^i film 
in the IN2 interface. 

7. Thickness t/L as a function of the undersaturation 
e = 1 — p/pt from the triple point pressure pt for 
diameter ratios d = 4.0(o), 4.2(o), 4.5(o), 5.0(>) of 
thin-thick mixtures. The inset shows the film thick- 
ness t/Li for long-short mixtures with length ratio 
I = 3.0(0), 3.1(0). 



8. Correlation length for rods in the triple-point Ni 
phase of thin-thick mixtures as a function of diam- 
eter ratio d, determined from adsorption (*) analy- 
sis (Eq. Ill|l and from density asymptotics (o) (Eq. 
I12|) . Inset (a) shows In |(5/3(z, 6')| for several values 
of 9 at the A'^i side of the INi interface for d — 4.0. 
Inset (b) displays the correlation length (in units of 
length of short rods Li) for long-short mixtures as 
a function of length ratio I. The dashed lines indi- 
cate the correlation length for rods in the N phase 
of the monodisperse hard-rod fluid. 

9. Surface tension ratio R [sec Eq. (|13|l ] as a function 
of the triple point under-saturation e for diame- 
ter ratio d = 4.0(o), 4.2(o), 4.5(<i), 5.0(i>). The inset 
shows the same quantity for long-short mixtures 
/ = 3.0(0), 3.1(0). 
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FIG. 1: K. Shundyak and R. van Roij 
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FIG. 4: K. Shundyak and R. van Roij 

FIG. 3: K. Shundyak and R. van Roij 
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FIG. 8: K. Shundyak and R. van Roij 



